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ABSTRACT 

A  gas  turbine  can  combustor  is  designed  to  bum  the  fuel  efficiently  while  reducing  the  NOx 
and  CO  emissions,  and  lowering  the  wall  temperature.  Environmental  challenges  with  gas  turbine 
include  low  levels  of  NOx,  CO  and  soot  amongst  other  pollutants.  There  is  a  need  for  new  concepts 
and  technology  to  satisfy  the  pollutants  emission  regulations  and  to  enhance  energy  conservation. 
Specifically,  ultra-low  NOx  combustor  technology  is  required  to  meet  the  ozone  depletion  challenge. 
Researchers  now  face  a  challenge  of  developing  dry  low-NOx  emitting  stationary  and  aero  engines. 
However,  any  concept  for  environmental  pollution  control  requires  a  detailed  understanding  of  the 
physical  and  the  chemical  processes  that  occur  during  combustion.  In  this  paper,  a  three  dimensional 
numerical  investigation  of  the  Combustion  methane  air  mixture  in  a  gas  turbine  can  combustor  is 
carried  out  by  using  ANSYS.  The  objective  of  the  study  is  to  understand  the  combustion  phenomena 
at  different  planes.  The  various  parameters  like  air-fuel  ratio,  velocity  of  primary  air  inlet  are  used  to 
investigate  the  effects  on  parameters  like  combustion  chamber  on  different  plane  performance  and 
emission.  A  premixing  tube  is  augmented  with  the  combustion  chamber  which  has  a  primary  air  inlet 
port  and  three  gaseous  fuel  inlet  ports.  Air-methane  mixture  is  considered  to  enter  the  combustion 
zone  with  inlet  swirl.  The  homogeneity  of  mixture  before  and  after  swirl  and  other  important 
conditions  are  calculated  from  simulation  and  are  reported  with  the  help  ANSYS  FLUENT.  Weighted 
averages  of  velocity  magnitude  distribution  and  mass  fractions  of  methane  have  been  studied  at 
different  planes. 

Keywords:  Air-methane  mixture,  ANSYS  FLUENT,  Combustion  chamber,  Gas  turbine,  Can 
combustor. 


1.  INTRODUCTION 
1.1.  Gas  Turbine 

A  Gas  Turbine  is  a  heat  engine  which 
uses  fuel  energy  to  produce  mechanical  output 
power,  either  as  torque  through  a  rotating  shaft 
(industrial  gas  turbines)  or  as  jet  power  in  the 
form  of  velocity  through  an  exhaust  nozzle 
(aircraft  jet  engines),  thus  producing  electrical 
energy  from  burning  a  combustible  mixture  of 
fuel  (e.g.  natural  gas  or  evaporated 
hydrocarbons)  and  air.  When  the  gas  mixture 
burns,  the  volume  of  the  gas  increases.  This 
expansion  in  gas  volume  makes  the  rotor  of  the 
turbine  to  rotate  and  this  rotation  is  then 


converted  into  electrical  energy.  The 
schematic  diagram  of  a  gas  turbine  is  shown  in 
figure  Al.  There  are  two  important  families  of 
gas  turbines: 

1.  Stationary  gas  turbines:  These  type 
of  turbines  are  used  to  produce 
power  in  large  scales,  for  example 
in  power  plants. 

2.  Turbofan  and  turbojet  gas 

turbines:  These  turbines  are 

usually  used  as  aircraft  engines 
and  are  sometimes  referred  to  as 
jet  engines.  A  variety  of  turbofan 
and  turbojet  gas  turbines  are  used 
in  military  and  commercial 
aircrafts. 
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1.2.  Combustion  chamber 

The  combustion  chamber  is  the  place 
where  two  major  events  take  place;  one  is  the 
mixing  of  air  and  fuel  and  the  second  is  the 
actual  combustion.  Depending  upon  whether 
premixed  flame  combustors.  A  brief 
explanation  about  the  two  types  of  combustors 
is  given  below. 

•  Diffusion  Flame  Combustors:  In 
diffusion  flame  combustors,  fuel  and  air 
mixing  and  combustion  takes  place 
simultaneously.  Speed  of  flame  is 
limited  by  the  rate  of  diffusion.  These 
kinds  of  combustors  are  simple  to  build 
and  operate,  but  they  are  not 
environmentally  clean.  The  major 
drawback  of  these  combustors  is  that 
the  flame  exists  mainly  at  stochiometric 
conditions.  This  can  result  in  high  rates 
of  NOx  production. 

•  Premixed  flame  combustors:  These 
combustors  are  newer  than  the  diffusion 
flame  combustors.  They  mix  the  fuel 
with  air  to  a  high  degree  thus  the  flame 
exists  where  the  fuel  exists,  if  it  can  be 
stabilized.  Contrary  to  the  previous 
combustor  type  these  combustors  are 
more  complex  and  harder  to  design,  but 
they  produce  less  NOx.  Every  day  new 
challenges  arise  for  gas  turbines. 
Different  factors  like  increase  in  oil 
price,  new  type  of  fuels  like  bio  fuels, 
different  designs  like  premix 
combustors  and  many  other  factors  will 
challenge  engineers  to  develop  new 
combustors  or  improve  the  existing 
combustors. 

However,  in  order  to  achieve  a  smooth 
burning,  air  and  fuel  should  be  mixed  before 
burning.  During  the  combustion  the  fuel-air 
mixture  ignites  and  raises  the  temperature. 
Rise  in  temperature  increases  the  volume 
which  in  turn  drives  the  fluid  forward.  The 
design  of  the  combustion  chamber  plays  a  vital 
role  in  controlling  the  combustion  process. 
Some  important  problems  occurring  as  a  result 
of  poor  design  of  a  combustion  chamber  are 
listed  below: 

1)  Poor  mixing:  When  fuel  is  not  mixed 
enough  with  air,  it  can  bum 
incompletely  which  results  in 
increased  levels  of  CO,  soot,  NOx  and 
unburned  hydrocarbons  (UHC) 

2)  Uneven  combustion:  This  happens 
when  temperature  of  a  section  goes 


the  mixing  of  air  and  fuel  is  done  before  the 
combustion  or  during  the  combustion, 
combustors  are  divided  into  two  groups 
namely  diffusion  flame  combustors  and 

high  while  the  neighboring  sections 
are  colder,  which  can  result  in  extra 
thermal  stresses.  At  times  thermal 
stresses  may  lead  to  material  fatigue 
and  failure. 

3)  Environment:  Incompletely  burned 
gases  or  unburned  hydrocarbons 
(UHC)  can  poison  the  environment. 
UHC,  NOx  and  soot  are  important 
factors  for  each  burning  device.  The 
design  should  lower  them  as  much  as 
possible. 

4)  Economy:  With  increasing  price  of  oil, 
it  is  important  that  gas  turbines  have 
high  efficiency  and  therefore  low  fuel 
consumption.  One  of  the  most 
important  parts,  in  order  to  achieve 
high  efficiency,  is  the  combustion 
chamber. 

1.3.  Previous  work  in  the  area 

The  development  of  advanced 
combustion  capability  for  gaseous  hydrogen 
and  hydrogen-blended  hydrocarbon  fuels  in 
gas  turbine  application  is  an  area  of  much 
interest.  In  most  countries  NOx  emissions  of 
stationary  gas  turbines  are  regulated  strictly. 
To  satisfy  these  regulations,  modem  gas 
turbine  manufacturers  rely  mostly  on  lean  pre¬ 
mixed  combustion.  Mixing  large  amounts  of 
air  with  gas  prior  to  burning  reduces  the  flame 
peak  temperature  and  thereby  leads  to  lower 
NOx  and  CO  emissions.  Efficient  combustion 
is  the  key  integral  part  not  only  in  power 
generation  related  sectors  like  gas  turbine, 
aircraft  engines  etc.  but  also  in  industrial 
burners  where  in  diffusion  flame  has  been  in 
use  for  a  long  time. 

In  the  present  study,  emphasis  is 
placed  on  the  fuel  preparation  process  which  is 
an  important  factor  and  sometimes  play  a  key 
role  in  the  reduction  of  pollutant  emission  from 
gas  turbines  [1].  It  is  a  challenge  to  handle  the 
environmental  problems,  with  low  pollution 
from  a  broad  range  of  fuels  [2].  New  concepts 
introduced  to  gas  turbine  industry  include  lean 
premixed  combustion  (LPM),  rich-burn  quick- 
Quench  lean-burn  (RQL)  combustion,  catalytic 
combustion  and  selective  catalytic  reduction 
(SQR);  high  -  efficiency  low  emissions 
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are  being  extended  to  non-premium  fuels 
such  as  coal  gas  and  crude  oil;  new  materials 
such  as  super  alloys,  thermal  barrier  coatings, 
and  ceramics  have  been  incorporated  into 
designs;  and  improved  theories  are  greatly 
dependent  on  advanced  laser-based  diagnostics 
of  flame  structure  and  have  led  to  design  tools 
of  increasing  scope  [3].  Among  these,  the  lean 
premised  combustion  appears  to  be  the  most 
effective  mainly  because  of  its  two -fold 
advantages  [4]  and  [5]  i.e.,  Reduction  of  flame 
temperature  owing  to  the  presence  of  excess 
air  which  eliminates  production  of  thermal 
NOx  and  low  exhaust  temperature  thereby 
increasing  the  life  time  of  turbine  blades  and 
other  mechanical  components. 

The  primary  drawback  of  such  LPM 
combustors  is  reduction  in  the  operability 
range  between  the  flashback  and  blowout 
regime  which  decreases  the  overall  stability  of 
the  combustion  process  [6].  Flame  stabilization 
schemes  that  have  been  widely  used  in  such 
applications  mainly  involve  bluff  body  or  swirl 
stabilization  or  a  combination  of  both.  The 
physical  process  that  anchors  a  premixed  flame 
behind  a  bluff -body  flame  holder  involves 
mixing  of  the  combustible  mixture,  which  are 
residing  in  the  recirculation  zone.  The  popular 
stabilization  mechanisms  used  for  sustaining 
the  flame  in  the  combustor  are  bluff  body 
flame  holders  which  is  reported  by  the 
previous  researchers  [7]  and  are  taken  into 
account  by  the  present  research  group.  In  this 
process  bluff-body  flame  holder  involves  in 
mixing  of  the  combustible  mixture  with  the  hot 
combustion  products  residing  in  the 
recirculation  zone,  thus  facilitating  continuous 
ignition  of  the  premixed  reactants  in  its 
vicinity  are  considered  and  used  in  the  process 
used  by  the  present  research  group  as  a 
parameter.  The  popular  stabilization 
mechanisms  used  for  sustaining  the  flame  in 
the  combustor  are  bluff  body  flame  holders  [8]. 
Dump  plane  swirler  and  pilot  flame  Lean 
blowout  (LBO)  is  a  very  important  aspect  of 
combustion  performance  for  gas  turbine 
combustor.  The  problem  of  LBO  was  always 
regarded  as  relatively  minor.  This  was  to  a 
large  extent,  due  to  the  almost  universal  use  of 
pressure-swirl  spray  atomizers.  The  poor 
mixing  characteristic  of  these  pressure  swirl 
atomizer  creates  many  performance  problems, 
the  last  being  a  high  rate  of  soot  formation  in 
the  primary  zone.  However,  this  poor  mixing 
of  fuel  and  primary  air  has  the  useful 


advantage  of  allowing  combustion  to  occur.  In 
recent  years,  the  trend  towards  improving 
fuel/air  mixture  prior  to  combustion  in  order  to 
reduce  pollutant  emission  has  led  to  a 
narrowing  of  stability  limits  and  increasing  the 
satisfactory  LBO  performance.  A  failure  in 
sustaining  the  flame  in  the  combustors  lead  to 
the  phenomenon  called  lean  blowout  (LBO). 
To  avoid  such  unwanted  phenomena,  current 
combustion  are  typically  operated  with  a  wide 
margin  above  the  lean  blowout  limit  which 
tends  to  increase  the  NOx  formation  from  the 
engine.  Hence,  a  significant  optimized  design, 
balancing  operation  economy  and  reduction  of 
emission  is  crucial.  The  current  work  is  aimed 
at  addressing  a  part  of  this  issue  where  a  CFD 
analysis  is  carried  out  on  dump  plane  swirler 
geometry  to  investigate  the  homogeneity  of 
air-fuel  mixture  [9].  A  complete  active  control 
system  sensing,  actuation,  and  control 
algorithm  has  been  developed  to  permit  turbine 
engine  combustion  to  operate  safely  closer  to 
the  lean-blowout  (LBO)  limit,  even  in  the 
presence  of  disturbances.  The  system  use  OH 
chemiluminescence.  Chemiluminescence  is 
directly  related  to  chemical  reaction  rates  and 
it  can  provide  information  on  the  presence  and 
strength  of  the  combustion  process  in  a 
specific  region  of  the  combustor  [10].  Though 
plenty  of  sensing  tragedies  for  predicting  the 
proximity  of  combustion  to  lean  blowout  have 
been  reported  in  the  past,  the  general  trend  is  to 
study  the  flame  dynamic  near  the  blowout  limit 
at  a  fixed  degree  of  premixedness.  The  degree 
of  premixedness  influence  heat  release 
fluctuations  which  being  the  primary  cause  of 
dynamic  instability  and  subsequently  affect  the 
LBO  limit  [11].  Several  experimental  studies 
characterizing  the  blowout  phenomenon  along 
with  secondary  measurement  like  strain  rate  at 
flame  base  pulse  rate  before  the  onset  of  lift  - 
off  have  been  done  using  modem  optical 
diagnostic  tools  like  LIPF(Laser-induced 
Predissociative  Fluorescence  ),  PIV(Particle- 
Image  Velocity)  and  OH  PLIF(Planar  Laser 
Induced  Fluorescence)  [12].  Optical  OH 
chemiluminescence  has  also  been  used  by  [13] 
to  observe  extinction  and  re-ignition  events 
near  the  blowout  event  which  has  been  called 
as  precursor  events.  LBO  prediction  is  a  very 
important  aspect  in  combustion  scenario  and 
various  methods  have  been  established  for  the 
same  over  the  years.  The  very  primitive  of 
such  methods  was  based  on  monitoring 
pressure  in  the  combustion  chamber. 
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1.4.  Scope  of  present  work 

The  focus  of  the  present  work  is  to 
study  the  dynamic  characteristics  of  fuel-air 
mixing  in  a  model  gas  turbine  dump 
combustor.  The  model  geometry  is  shown  in 
Fig.  2,  where  air  is  entering  at  one  end  of  the 
combustor  and  gaseous  fuel  is  entering  from 
one  of  the  perpendicular  ports  shown  (port-1). 
The  air  velocity  is  varied  sinusoidally  with 
time  and  the  effect  of  this  mixing  parameter  is 
studied  at  different  planes.  The  work  is  carried 
out  numerically  using  FLUENT  software. 

2.  CALCULATION  METHODOLOGY 

2.1.  Problem  description 

The  combustion  chamber,  3D 
geometry,  is  generated  using  ANSYS 
FLUENT.  The  combustion  zone  is  cylindrical 
in  shape  with  20cm  length  and  6cm  internal 
diameter.  Before  the  combustion  zone  a 
premixing  tube  is  provided  with  a  maximum 
premixing  length  of  35cm  and  diameter  of  2.3 
cm.  Dump  plane  is  considered  as  the  plane  at 
which  the  combustion  zone  starts  with  an 
abrupt  increase  of  cross-section.  At  the  bottom 
of  the  premixing  tube  i.e.  37.5  cm  from  the 
dump  plane  the  air  inlet  port  of  1.5  cm  diameter 
is  provided.  There  are  3  fuel  inlet  ports  each  of 
6mm  diameter  at  10  cm  distance  from  each 
other,  the  uppermost  being  15cm  below  the 
dump  plane.  The  fuel  inlet  ports  which  are  at  a 
distance  of  35cm,  25cm  and  15cm  from  the 
dump  plane  are  named  as  port  1,  port  2  and 
port  3  respectively.  Just  below  the  dump  plane 
a  swirler  of  length  1.5  cm  and  blade  angle  with 
axial  direction  60°  is  placed.  The  number  of 
blades  on  the  swirler  is  8  and  blade  thickness  is 
1mm.  Diameter  of  the  swirler  section  is  2.9cm. 
Mesh  generated  is  unstructured  and  in  standard 
form  as  provided  by  ANSYS.  The  number  of 
nodes  and  elements  are  337  land  12346 
respectively.  Details  may  be  found  in  figure 
A2. 

For  Numerical  simulation  of  the 
generated  meshed  geometry  in  3D,  pressure- 
based,  transient,  species  transport  and  k-e 
turbulent  models  are  utilized.  For  species 
transport  model  methane-air  mixture  is  used 
for  all  cases.  For  fuel  (Methane),  entry  through 
port-1  the  air  flow  rate  of  85  litres  per  minute 
(LPM)  are  used,  the  corresponding  air  inlet 
velocity  being  8.0208m/s.  For  entry  of  fuel 
(Methane)  through  port  -1,  the  premixing 
length  is  35  cm.  The  fuel  flow  rate  is  varied 


and  equivalent  ratio  and  fuel  entry  port  is 
considered  as  1  for  each  case  respectively.  For 
this  equivalent  ratio  the  fuel  inlet  velocity  is 
maintained  at  1.062m/s.  The  most  common 
combustors  have  separate  inlets  for  fuel  and 
air.  Both  the  fuel  and  the  air  are  assumed  to  be 
perfectly  mixed  before  combustion  chamber. 
Homogeneity  in  air-fuel  mixture  with  varying 
premixing  length  with  respect  to  contours  of 
mass  fraction  of  methane  (CH4)  and  velocity 
magnitude  at  different  important  planes  inside 
and  outside  the  combustor  is  found.  Planes  are 
generated  at  a  distance  of  3.5cm  before  the 
dump  plane  and  15cm,  and  20  cm  after  dump 
plane.FigureA3  and  figureA4  show  the  graphs 
and  contour  constructed  to  display  the 
distribution  of  mass  fraction  of  methane  and 
velocity  magnitude  along  a  time.  Results  are 
presented  for  most  important  comparisons  and 
variations  with  time. 

2.2.  Boundary  conditions 

Boundary  conditions  are  as  follows 

1)  Air  Inlet  port:  UDF  Unsteady 
Velocity, 

2)  Turbulence:  Specification  Method  - 
Intensity  and  Hydraulic  Diameter. 
Turbulent  intensity  -10% 

Hydraulic  Diameter-0.0 15m 

3)  Port-1  :  Velocity  1.062(m/s) 
Turbulence:  Specification  Method  - 
Intensity  and 

Hydraulic  Diameter.  Turbulent 
intensity  -10% 

Hydraulic  Diameter-0.06m 

4)  Walls:  Adiabatic  and  free  slip 

5)  Reference  pressure:  0  (Pascal) 

2.3.  Methodology 

A  commercial  available  CFD  code 
FLUENT  has  been  used  for  the  analysis.  Five 
equations  are  used  to  model  the  flow  field. 
These  equations  are  continuity,  momentum, 
swirler  velocity,  species  transport  and 
turbulence.  The  following  sections  will 
describe  these  equations. 

2.4.  UDF  equation 

The  UDF  equation  is  shown  in  equation  (2.1) 

^  +  V.(pv)  =  0  (2.1) 

2.5.  Continuity  equation  or  mass 
conservation  equation 
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This  equation  describes  that  what 
comes  in  should  go  out.  Mass  can  neither  be 
destroyed  nor  created. 

The  general  form  of  the  mass 
conservation  equation  is  given  by  the  e 
following  equation 


^  +  V.  (pvv)  -  -Vp  +  V.  (x  )  +  pg  +  F  (2.2) 
Equation  (2.2)  is  and  is  valid  for 
incompressible  as  well  as  compressible  flows. 


2.6.  Momentum  conservation  equation 

Momentum  is  a  vector  quantity  that  is 
the  product  of  mass  by  velocity  vector.  In  a 
closed  system,  momentum  cannot  be  created 
nor  destroyed.  It  should  be  conserved.  The 
equation  (2.3)  reads, 

f  =  p[(Vv  +  VvT)-\  V.vl]  (2.3) 


Where  p  is  the  static  pressure,  r  is  the 
stress  tensor  (described  below),  mdpg  and  F 
are  the  gravitational  body  force  and  external 
body  forces,  respectively.  F  also  contains 
other  model-dependent  source  terms  such  as 
user-defined  sources.  The  stress  tensor  t  is 
given  by  equation  (2.4). 


yt(pw)  +  1-fx(rpuw)  +  1-  £  ( rpvw )  = 
[rndf}  +  4 f[rVf  (-)]-—  ■  (2.4) 

r  dx L  dxJ  r  9rL  drKryj  r 


determined  as  one  minus  the  sum  of  the  N  -  1 
solved  mass  fractions.  To  minimize  numerical 
error,  the  Nth  species  should  be  selected  as  that 
species  with  the  overall  largest  mass  fraction, 
such  as  N2  when  the  oxidizer  is  air. 

2.9.  Turbulence  models 

2.9.1.  Mass  diffusion  in  turbulent  flows 

In  turbulent  flows,  ANSYS  FLUENT 
computes  the  mass  diffusion  equation  (2.7)  in 
the  following  form, 

A  =  -  (  pDi,m  +  )  Vy;  -  DT,i  Y  (2.7) 

Where  Set  is  the  turbulent  Schmidt  number,  pt 
is  the  turbulent  viscosity  and  Dt  is  the 
turbulent  diffusivity.  The  default  is  Set  0.7. 

2.9.2.  k -e  Turbulence  model 

The  k-£  model  is  one  of  the  most 
common  turbulence  models.  It  is  a  two 
equation  model  which  means  that  two  extra 
transport  equation  is  included  to  represent  the 
turbulent  properties  of  the  flow.  In  these  cases, 
using  mixture  properties  and  mixture  velocities 
is  sufficient  to  capture  important  features  of 
the  turbulent  flow.  The  k  and  £  equations  (2.8) 
and  (2.9)  describing  this  model  are  as  follows 
^(pmk)  +  V.  (pmv^k)  =  V.  Vk)  +  Gkm  -  pme  (2.8) 


2.7.  Momentum  conservation  equation  for 
swirl  velocity 

The  tangential  momentum  equation 
(2.5)  for  2D  swirling  flows  may  be  written  as 

£  ( PVi)  +  V.  ipvyt)  =  -V.h  +  Ri+St  (2.5) 


2.8.  Species  transport  equations 

When  you  choose  to  solve 
conservation  equations  for  chemical  species, 
ANSYS  FLUENT  predicts  the  local  mass 
fraction  of  each  species,  Yi,  through  the 
solution  of  a  convection  diffusion 
equation(2.6)  for  the  ith  species. 

^  (pyO  +  V.  (pvyO  =  -V.fl  +  Ri+  Si  (2.6) 

Where  Ri  is  the  net  rate  of  production  of 
species,  i  be  the  chemical  reaction  and  Si  is  the 
rate  of  creation  by  addition  from  the  dispersed 
phase  plus  any  user-defined  sources.  In  this 
equation  (2.6)  we  assume  that,  Si  is  equal  to 
zero.  An  equation  of  this  form  will  be  solved 
for  N  -1  species  where  N  is  the  total  number 
of  fluid  phase  chemical  species  present  in  the 
system.  Since  the  mass  fraction  of  the  species 
must  sum  to  unity,  the  Nth  mass  fraction  is 


jt  (pme)  +  V.  (pmv^e)  =  V.  (  ^  Ve  )  + 
l(CleGKm-C2€pm6 )  (2.9) 

where  the  mixture  density  and  velocity,  pm and 
vm  are  computed  from  eqautions(2.10)  and 
(2.11) 

Pm~^i=i^iPi  (2.10) 

and 


- *  ^i=iaiPiv: 

Vm  =  Sf=1«iPi 


(2.11) 


The  turbulent  viscosity,  m  is  computed  from 
equation  (2.12) 

Pt,m  ~  PmCp~  (2.12) 

and  the  production  of  turbulence  kinetic 
energy,  Gk  m ,  is  computed  from  equation(2.13) 


Gk,m  =  (V%)T}:V%  (2.13) 


3.  RESULTS  AND  DISCUSSIONS 

The  area  weighted  average  of  methane, 
mass  fraction  of  CH4,  area  weighted  average 
velocity  magnitude,  mass  weighted  average  of 
methane  are  calculated.  Figure  5  shows  the 
Mass  weighted  Average  Velocity  magnitude 
contour  plots,  where  it  is  found  that  the  result 
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for  dump  plane,  dump  plane  attached  with 
s wirier  and  the  plane  created  are  3.5cm,  15  cm 
and  20  cm  respectively.  For  air  flow,  velocity 
of  UDF  [where  uo  =8.062,  A  =  4%  of  uo  w  = 
95]  at  air  inlet  port  (port  -1)  having  a  constant 
velocity  is  1.062m/s.  While  Fig.  8  shows  mass 
weighted  average  velocity  magnitude 
distribution  along  the  above  mentioned  planes 
with  respect  to  time. 

The  scale  utilized  in  important  planes 
before  and  after  dump  plane  is  0  to  5m/s  which 
is  significant  with  respect  to  the  Velocity 
Magnitude.  For  the  planes  constructed  before 
and  after  dump  plane  the  same  scale  is  used 
and  methane  mass  fraction  of  CH4  scales  0  to 
0.023. 

3.1.  Velocity  magnitude 

Homogeneity  in  air-fuel  mixture  is 
obtained  by  using  the  UDF  Velocity  at  air  inlet 
port  and  constant  velocity  at  port- 1. The 
variation  of  velocity  magnitude  with  respect  to 
distance  and  time  in  different  important  plane 
before  swirler  and  after  swirler  is  obtained.  We 
study  from  the  contour  at  YZ  plane  (figureA5) 
and  graph  (figureA7)  that  the  velocity 
magnitude  is  very  large  before  swirler  and 
after  swirler.  Attached  with  the  swirler  is 
called  dump  plane.  Another  two  important 
plane  in  a  combustion  chamber  is  15  cm  after 
dump  plane  and  20cm  after  dump  plane  in 
which  the  velocity  decreases  gradually  with 
distance  as  well  as  time  as  the  flow  progresses 
through  the  swirler(magnitude  of  velocity 
before  swirler  plane  is  3.5  m/sec.  and  after 
swirler  dump  plane  it  is  0.56  m/sec„15cm  after 
dump  plane  it  is  0.494  m/sec„20cm  after  dump 
plane  it  is  0.479  m/sec)  .After  swirler,  flow 
enters  into  the  combustion  chamber  inlet  that  is 
called  the  dump  plane  where  velocity  further 
reduces  due  to  large  diameter  difference 
between  before  swirler  plane  (<p2.3cm)and 
after  swirler  plane(<p6cm)  resulting  in  the 
formation  of  wall  recirculation  zone 
(figure  A6).  The  plane  diameter  in  a 
combustion  zone  is  more  than  double  with 
respect  to  premixed  tube  diameter  and  swirlers 
also  play  a  measure  role  to  oppose  the  velocity 
of  Air-fuel  mixture. 

3.2.  CH4  mass  fraction 

In  order  to  achieve  better  mixing 
between  fuel  and  air  in  combustor,  turbulent 
flow  must  be  generated  to  promote  mixing. 
Below  figure  shows  the  CH4  mass  fraction 
contours  at  YZ  planes  while  figureA9(a), 


shows  CH4  mass  fraction  contours  at  plane 
before  swirler  compared  with  dump  plane  and 
other  two  different  important  planes  after 
swirler  .FigureA9  is  important  as  they  show 
homogeneity  of  air  -fuel  mixture  and  mixture 
pattern  at  different  position  along  the  length  of 
the  combustor.  Comparing  (a),  (b),  (c),  (d)  of 
figureA9,  it  can  be  easily  detected  that  in  the 
mixture  before  swirler  and  after  swirler  in  the 
combustor  the  methane  mass  fraction  regularly 
decreases  at  different  plane  with  increase  in 
distance  and  time.  FigureAlO  shows  a  plot  of 
mass  fraction  of  CH4  versus  time  for  different 
planes.  There  was  decreasing  fluctuation  of 
mass  fraction  of  CH4  while  increasing  the 
distance.  The  fluctuation  difference  between 
maximum  and  minimum  methane  mass 
fraction  values  are  0.0006,  0.0001,  0.000022, 
0.0000038  for  plane  position  3.5cm  before 
dump  plane,  dump  plane,  15cm  after  dump 
plane,  20cm  after  dump  plane  respectively. 
This  shows  more  rapid  decrease  of  the 
difference  for  plane  20cm  after  dump  plane 
when  compared  to  3.5  cm  before  dump  plane. 
This  means  that  fluctuation  is  almost  steady  on 
the  plane  20cm  after  dump  plane  with  respect 
to  3.5cm  before  dump  plane.  As  the  premixed 
length  increases  homogeneity  enhances  in  the 
mixture.  With  variation  of  distance  the  effect 
on  homogeneity  is  dependent  on  the  premixed 
length.  FigureAll  shows  the  Mass-weighted 
Average  Mass  fraction  of  CH4  distribution 
along  different  planes  with  respect  to  time 
while  figureA12  to  16  depict  the  Area- 
Weighted  Average  Velocity  magnitude 
distribution  for  different  amplitudes  for  W-  5, 
W=  20,  W=35,  W=50,  W=80  and  W=  95 
respectively. 

4.  CONCLUSIONS 

The  gas  turbine  combustor 
development  will  have  greater  focus  on  both 
environmental  pollution  and  energy 
conservation.  Advanced  concepts  for 
environment  pollution  control  required  detailed 
understanding  of  the  physical  and  chemical 
process  occuring  during  combustion.  Test  data 
on  laboratory  Scale,  bench  scale,  and  prototype 
combustor  are  urgently  needed  for  determining 
the  details  of  the  chemistry  involved,  fuel-air 
mixture,  fuel  injector  design,  combustor 
geometry,  combustion  air  swirl,  flow 
distribution  and  fuel  spray  characteristics.  It 
should  be  mentioned  in  the  velocity  and  mass 
fraction  of  CH4  plots  that  mass  fraction  of  CH4 
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and  velocity  distribution  may  vary  with 
number  of  iteration  and  that  is  the  reason  for 
the  variastion  of  plots  and  figures.  The 
variation  of  the  mass  fraction  of  methane  and 
velocity  fall  into  acceptable  error  range  for 
numerical  error,  thus  the  results  were  accepted. 
The  numerical  simulation  and  the  results 
qualitatively  determine  how  the  homogeneity 
of  air-fuel  mixture  in  a  gas  turbine  combustor 
varied  with  variation  of  distance  as  well  as 
time.  Dump  plane  results  are  extremely 
important  in  this  aspect  as  it  signifies  the  first 
plane  after  the  mixtures  enter  the  main 
combustion  zone.  As  the  premixing  length 
increases  homogeneity  enhances  in  the 
mixture.  With  variation  of  distance  the  effect 
on  homogeneity  depends  on  the  premixing 
length. 

5.  Scope  of  future  work 

The  dynamic  characterization  of 
mixing  can  be  studied  with  different  wave 
forms.  Dynamic  characteristic  of  combustion 
can  also  be  studied.  System  identification  foe 
mixing  as  well  as  combustion  can  be  carried 
out. 
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APPENDIX  A 
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Figure  Al. Schematic  of  gas  turbine  (P-1) 
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Figure  A2. Meshed  geometry  used  for  ANSYS  FLUENT  simulation 
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Figure  A3. Contour  of  mass  fraction  of  methane  (CH4)  at  plane  created  before  and  after  dump  plane  by  using 

ANSYS  FLUENT 
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Figure  A4. Contour  of  velocity  magnitude  at  Plane  created  before  and  after  dump  plane  by  using  ANSYS 

FLUENT 
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Figure  A5.  Velocity  Magnitude  Contour  on  Different  plane  for  entry  of  fuel  through  port-1  for  velocity 
V=1.062m/s  and  air  flow  Velocity  of  UDF  (a)3.5  cm  before  dump  plane  (b)dump  plane  (c)  15cm  after  dump 

plane  (d)  20cm  after  dump  plane. 


Figure  A6. Recirculation  zone  and  s wirier 
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Figure  A7. Area- weighted  average  velocity  magnitude  distribution  along  different  planes  with  respect  to  time 
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Figure  A8. Mass-weighted  average  velocity  magnitude  distribution  along  different  planes  with  respect  to  time 
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Figure  A9.Mass  fraction  of  CH4  contour  on  different  plane  for  entry  of  fuel  through  port-1  for  velocity 
V=1.062m/s  and  air  flow  Velocity  of  UDF  (a)3.5  cm  before  dump  plane  (b)dump  plane  (c)  15cm  after  dump 
plane  (d)  20cm  after  dump  plane 
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Figure  A10.  Area- weighted  average  mass  fraction  of  CH4  distribution  along  different  planes  with  respect  to  time 
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Figure  A1 1. Mass-weighted  Average  Mass  fraction  of  CH4  distribution  along  different  planes  with  respect  to  time 
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Figure  A1 2.  Area- weighted  average  velocity  magnitude  distribution  for  different  amplitudes  where  W  =  5 
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Figure  A 13. Area- weighted  average  velocity  magnitude  distribution  for  different  amplitudes  where  W  =  20 
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Figure  A1 4.  Area- weighted  average  velocity  magnitude  distribution  for  different  amplitudes  where  W  = 
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Figure  A 15. Area- weighted  average  velocity  magnitude  distribution  for  different  amplitudes  where  W  =  50 


Figure  A1 6.  Area- weighted  average  velocity  magnitude  distribution  for  different  amplitudes  where  W  =  80 
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